Polarization contributions to the spin-dependence of the effective interaction in 

neutron matter 
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We calculate the modification of the effective interaction of particles on the Fermi surface due 
to polarization contributions, with particular attention to spin-dependent forces. In addition to 
the standard spin-spin, tensor and spin-orbit forces, spin non-conserving effective interactions are 
induced by screening in the particle-hole channels. Furthermore, a novel long-wavelength tensor 
force is generated. We compute the polarization contributions to second order in the low-momentum 
interaction Vi ov ,k and find that the medium-induced spin-orbit interaction leads to a reduction of 
the 3 P2 pairing gap for neutrons in the interior of neutron stars. 



Introduction. - Landau-Fermi liquid theory is a pow- 
erful effective theory for strongly interacting Fermi sys- 
tems at low temperatures. It has been successfully ap- 
plied to liquid 3 He, nuclear matter and nuclei. While 
the free interaction between 3 He atoms is almost state- 
independent, the nuclear interaction is complicated due 
to large non-central spin-orbit and tensor forces, which 
are crucial for understanding nuclear phenomena. For in- 
vestigations of matter under extreme conditions, such as 
nuclei with large proton or neutron excess and asymmet- 
ric nuclear matter in neutron stars, the role of non-central 
forces in the effective interaction must be understood. 

As part of a program to determine effective nuclear in- 
teractions using renormalization group methods yj, we 
analyze the spin-dependence of the quasiparticle inter- 
action and the low-energy scattering amplitude in the 
presence of non-central forces. We focus on pure neutron 
matter, and as an application, we estime the modifica- 
tion of the 3 P2 pairing gap in neutron star interiors due 
to the screening of the nucleon-nucleon interaction. This 
is a long-standing problem in neutron star structure, and 
since polarization effects suppress the S-wave gaps by a 
factor four large effects may be expected. 

Symmetry considerations. - In general, the two-body 
interaction is hermitian and constrained by symmetries, 
specifically, time-reversal and parity invariance, as well 
as invariance under exchange of particle labels. In ad- 
dition, in a non-relativistic theory in vacuum, the po- 
tential is Galilean invariant, i.e., independent of the 
particle-pair momentum P = pi + p 2 = P3 + P4- The 
possible operators are well-known |2j; scalar: 1, spin- 
spin: <Ti ■ er 2 , spin-orbit: i(<xi + <r 2 ) ■ q x q', tensor: 
S12 (q) = <r\ ■ qer 2 • q — 1/3 q 2 eri • er 2 (and the exchange 
thereof S\% (q')) as well as the quadratic spin-orbit force: 
qq' (ct 1 qcr 2 ■ q' + cri • q' er 2 • q- 2/3 q- q' <J\ -er 2 ). The 
momentum transfer is q = pi — P3 and in the exchange 
term q' = pi — p 4 . 

In the many-body medium the presence of the Fermi 
sea defines a preferred frame. Therefore, the effective 
two-body interaction depends on the two-body center of 
mass (cm) momentum. This is physically clear, since 



the magnitude of the cm momentum, for given momen- 
tum transfers, determines where the interacting particles 
are relative to the Fermi sea. For particles on the Fermi 
surface q 2 + q' 2 + P 2 = 4&§ and the momentum depen- 
dence of possible invariants is constrained geometrically 
since q, q' and P are orthogonal. As a consequence, the 
quadratic spin-orbit force vanishes in this case. 

Effective interactions on the Fermi surface. - In the 
presence of a Fermi sea, additional operators are possible. 
For particles on the Fermi surface, these are 

£12 (P) cm tensor (1) 

D 12 (q,P) = i(<ri - <x 2 ) q x P diff vector (2) 
j4i 2 (q',P) = (cr 1 x er 2 ) • (q' x P) cross vector (3) 

Operators Di 2 and An are related by exchange. These 
are antisymmetric in spin and thus do not conserve the 
spin of the interacting particle pair. In this Letter, we 
explore the microscopic origin of these forces and com- 
pute the contributions to the quasiparticle interaction in 
neutron matter. Our results can be used as input for 
calculations of neutron star properties. For particles not 
on the Fermi surface, further invariants are possible Q. 

Both the cm tensor and the cross vector operator sur- 
vive in the long-wavelength limit, q — > 0, and thus con- 
tribute to the quasiparticle interaction in nuclear matter. 
The exchange tensor in the quasiparticle interaction was 
considered in p, 6], and Landau parameters were com- 
puted in 0, • We introduce the Fermi liquid parame- 
ters, Hi, Ki, and Li for the non-central interactions, 
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' Pi =ff(cos0 q OSi 2 (qO + tf(cos0 q OSi2(P) 
P). 
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(4) 



The tensor operators are defined with unit vectors [2j 
and the dependence on Landau angle # q ' is expanded 
in Legendre polynomials, H(x) = J^i Hi Pi{x) etc. The 
novel Fermi-liquid interactions K and L have not been 
considered in previous work. Analyticity of the quasi- 
particle interaction implies H ~ q' 2 as q' — > 0, K ~ P 2 
as q' — * 2/cp (P — * 0) and L ~ q' as q' — * (as well as 
L ~ P as q' — > 2fep). Many-body effects may give rise 
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to singularities in the effective interaction, which modify 
these limits. This is illustrated by K (cos 6 4) in Fig. ^ 
where the pairing singularity cancels the zero of the cm 
tensor in the limit q' — > 2kp (cos 6^ — ► — 1). 

Effects of the many-body medium. - At second-order, 
there are contributions to the effective four-point vertex 
from scattering in the BCS channel with intermediate 
particle-particle and hole-hole excitations of cm momen- 
tum P, as well as scattering in the direct particle-hole 
or zero sound (ZS) channel and in the exchange particle- 
hole (ZS') channel. The ZS and ZS' channels include 
intermediate states with particle-hole excitations of mo- 
mentum q and q' respectively. For the quasiparticle in- 
teraction only the BCS and ZS' channel contribute. Gen- 
eral recoupling arguments imply that the interference of 
the spin-spin with the tensor force in the particle-hole 
channels leads to a large renormalization of the tensor 
interaction Q]. Moreover, the presence of a third particle 
in intermediate states induces novel contributions to the 
effective interaction, of the form of Eqs. © " © ■ 

For an antisymmetrized interaction / d (q, q') = 



, (q, q') - P< T V rTl ^ 2 (q', q) (where P„ is the spin ex- 



change operator and the superscript d labels the driving 
term pj), the particle- hole contributions in the ZS' chan- 
nel are given by 
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TABLE I: The I ^ 4 Fermi liquid parameters for neutron 
matter at the Fermi momentum fcp = 1.7 fin -1 , with effective 
mass m*/m = 0.83. The different contributions are discussed 
in the text. The total includes also the (small) boost correc- 
tions. The renormalization of the quasiparticle strength Zk F 
is neglected and we use Zk F = 1 throughout. 
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Landau / 





1 


2 


3 


4 


Vlow fc 


scalar F d 


-0.734 


-0.498 


-0.200 


-0.068 


-0.052 


spin-spin Gf 


0.842 


0.412 


0.219 


0.109 


0.053 


exch. tensor Hf 


0.529 


0.150 


-0.096 


-0.141 


-0.124 


ZS' channel 


scalar 


0.552 


0.406 


0.119 


0.131 


0.099 


spin-spin 


0.024 


-0.038 


-0.052 


-0.016 


0.002 


exch. tensor 


-0.214 


-0.218 


-0.086 


0.004 


0.067 


cm tensor 


-0.071 


-0.014 


0.104 


0.047 


0.009 


cross vector 


-0.015 


-0.073 


-0.057 


0.018 


0.032 


BCS channel 


scalar 


-0.291 


0.187 


-0.146 


0.121 


-0.088 


spin-spin 


0.032 


0.126 


-0.006 


0.032 


-0.019 


exch. tensor 


-0.201 


0.249 


-0.089 


0.073 


-0.060 


cm tensor 


-0.187 


0.160 


0.020 


0.002 


0.005 


total 


exch. tensor Hi 


0.077 


0.140 


-0.266 


-0.047 


-0.101 


cm tensor Ki 


-0.258 


0.147 


0.124 


0.048 


0.014 


cross vector Li 


-0.061 


-0.089 


-0.035 


0.025 


0.041 



where n p denotes the Fermi-Dirac distribution function 
and £ p is the quasiparticle energy for which a single- 
particle spectrum with effective mass is used. The ex- 
change tensor, iterated in the particle-hole ladder, yields 
an unusual ordering of the spin operators <T\ ■ t' <j\ ■ t <r 2 • 
ter 2 • t', where tW = ±(P/2 ± q/2 - p). This order- 
ing gives rise to a particular coupling between the spin 
and the angular motion, which leads to antisymmetric 
spin operators. Furthermore, particle-hole polarization 
contributions involving spin-orbit forces always result in 
spin non-conserving interactions. 

On the other hand, in the BCS channel the depen- 
dence on the cm momentum enters only through the 
phase space and the ordering of the spin operators is 
the same as in vacuum. Consequently, the BCS channel 
does not yield antisymmetric spin operators. The second 
order contribution is given by 
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where k,,/ = (q' ±q)/2 denotes initial/final relative mo- 
menta and /i is the chemical potential. 

Boost corrections. - In addition to the dynamical ef- 
fects, there are kinematical contributions from boosting 



the two-body interaction to the rest frame of the Fermi 
sea. To leading order in k^/m 2 , the boost to the frame 
where the nucleon pair carries momentum P is give n in 
terms of the vacuum interaction in the cm frame |l(i lll| 

^ CTl , CT2 (q,q',P) = (q,q') 

(<T 1 -<7 2 )xP-(q'-q)^ iff ,(q 1 q') 

2 Tt^M) -<ra) x P • (q' + q), (7) 



16 TO 2 

i 

16 771 



where T d = m* : k F f d j 'ir 2 . In Eq. J7J) the particles are 
restricted to the Fermi surface and direct and exchange 
terms are included. In the long-wavelength limit, we find 

<5^ CTl , CT2 (q',P) = --^ {F d (cos6V) 
4 vn 

+ G d (cos q » ) a x ■ <r 2 + H d (cos 9^ ) S 12 (q')} 

- -Ij Ai 2 (q', P){G d (cos# q + \ H d (cos# q 0}, (8) 
4m z D 

with standard notation for the scalar and spin-spin parts, 
F d and G d . We note that the boost corrections con- 
tribute to the spin non-conserving part, but not to the 
cm tensor Eq. ©. These kinematical effects are of or- 
der kp/m 2 , while the many-body effects are of order 
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FIG. 2: Comparison of the particle-hole induced spin-orbit 
amplitude to the free-space interaction. Results are presented 
for back-to-back scattering, P = 0, — 1.7 fm -1 (thick lines) 
and kp — 1.0 fm -1 (thin lines). Here, the (antisymmetrized) 
spin-orbit force is shown in units of the density of states with 
the operator Vso(cos # q ) i(ai + (T2) ■ q/&F x q'/kp. 



Fig. [21 wc find a repulsive induced spin-orbit interaction 
due to particle-hole screening. In order to qualitatively 
understand the resulting interaction, we assume a con- 
tact spin-spin and an averaged spin-orbit matrix element. 
Then, the ZS and ZS' channel for P = are repulsive 



FIG. 1: Angular dependence (cos# q / = 1 — q' 2 /2k F ) of the 
non-central quasiparticle interactions in neutron matter for 
&f = 1.7 fin . The Viowfc contribution as well as the polar- 
ization effects from the exchange particle-hole (ZS') and the 
particle-particle (BCS) channel are shown for m*/m = 0.83. 



(m*/m) fep(V). Here, the brackets denote the angular 
average over the interaction in the loop integral. 

Results. - We start from the free-space low-momentum 
interaction Viowfc E3 with a density-dependent cutoff 
A = V2k F 0. Detailed results for both neutron and nu- 
clear matter will be presented elsewhere. Here, we focus 
on the novel spin-dependent interactions. The contribu- 
tions to the Fermi liquid parameters are given in Tabled 
and the dependence on cos 6> q / is shown in Fig. ^ We 
find a substantial renormalization of the exchange ten- 
sor, which necessitates a self-consistent treatment within, 
e.g., the RG approach 0, and significant contributions to 
the new interaction terms. In particular, the cm tensor 
is comparable to the exchange tensor at this order. 

In calculations of transport processes and pairing phe- 
nomena, one needs the scattering amplitude at finite mo- 
mentum, q =/= 0. In particular, the in-medium modifica- 
tion of the spin-orbit interaction is of special interest, 
because it is crucial for an accurate reproduction of the 
P-waves in free space, and similarly for a realistic assess- 
ment of P-wave pairing in neutron stars. As shown in 



K£ d (cos0 q ) = -G d a — (V s d ) (U(q/k F ) + U(q'/k F )), 

(9) 

for an attractive spin-orbit force. Here U(q/k F ) de- 
notes the (positive) static Lindhard function and = 
k F =F 2 cos q . As in Fig. V^q is only weakly de- 
pendent on cos# q . The contributions due to the mixing 
of spin-orbit and tensor forces are also repulsive, with a 
similar but more complicated momentum dependence. 

Triplet pairing. - To illustrate the importance of 
non-central induced interactions, we estimate the angle- 
averaged gap using weak coupling BCS theory 13] (for 
details see The coupling to the 3 F2 partial wave 

is neglected in this exploratory calculation. In the weak 
coupling approximation, the 3 P2 gap is then given by 



A ^ F 

A3p 2 = — exp 



2fc F m^ owfc .3p 2 (fc F ,fc F ) 



(10) 



where the arguments of VWfc are the magnitude of 
the relative momenta |kj,/| = k F . When particle- hole 
screening effects are included, the pairing interaction 
V| owfc .3 P (k F ,k F ) is replaced by the 3 P2 projection of the 
effective interaction in the particle-particle channel. The 
high-lying states in the BCS channel are included in Viowfc 
while the low-lying part is accounted for through the (ap- 
proximate) solution of the gap equation. 

We find that spin-dependent polarization effects re- 
duce the 3 P2 pairing gap (C compared to A and B in 
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FIG. 3: The angle- averaged pairing gap A3p 2 in the 3 P2 chan- 
nel versus the Fermi momentum in neutron matter. The di- 
rect (V| ow fc) pairing gap, computed with the free and with an 
effective neutron mass, as well as the gap including particle- 
hole polarization effects on the pairing interaction are shown. 
We also give the gap, obtained when only the central and only 
the spin-orbit polarization contributions are taken into ac- 
count. For reference, we show the results of Baldo et al. flBf . 
obtained by solving the BCS gap equation in the coupled 
3 P2- 3 F2 channel for different free-space interactions. VIowfc is 
obtained from the CD Bonn potential. 

Fig. |3J). This is in contrast to the increase of the 3 P2 
gap, which one obtains when polarization effects due only 
to central forces are included (dotted vs. dashed line in 
Fig. |3J) [16| . The reduction of the gap is predominantly 
due to the repulsion from the medium-induced spin-orbit 
force, discussed above. This effect was not taken into ac- 
count in earlier work. Note that for the densities given in 
Fig- El the ratio of the induced pairing matrix element to 
the free-space Vj ow fc contribution ranges from 0.15 — 0.5. 
The significant reduction of the gap for only moderate 
changes of the pairing interaction is due to the singular 
dependence on the matrix element in Eq. l|10fl for small 
gaps. We emphasize that our results are qualitative; a 
quantitative calculation should include the full polariza- 
tion contributions for non-central interactions and the 
solution of the full coupled-channel BCS gap equation. 
We present results only for fcp < 2fm~ 1 , where the NN 
interaction is strongly constrained by data. 

Conclusions. - In summary, we have found novel 
non-central effective nuclear interactions, with spin non- 
conserving forces induced by particle-hole polarization 
effects. In microscopic shell model calculations, these 
are implicitly included in the polarization force of Kuo 



and Brown 17]. Furthermore, the renormalization of the 
spin-orbit interaction in the medium has important con- 
sequences for P-wave pairing. The resulting suppression 
of the superfluid gap has direct impact on the proper- 
ties of neutron stars and on their cooling by neutrino 
emission 0, 0] . The implications of the new interac- 
tions for nuclear spectra, for the spin-isospin response 
and neutrino transport in stellar collapse, for magnetic 
susceptibilities (see for deformations of the Fermi 
surface in spin-polarized systems, and for spin relaxation 
and mixing of spin and density waves remain to be in- 
vestigated. Since the new interactions contribute only 
to spin non-conserving transitions, it may be possible to 
observe these effects in scattering with polarized beams. 
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